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Abstract 

We demonstrate, by means of numerical simulation, a new mechanism to 
modulate and switch semiconductor lasers at THz and sub- r Hz frequency 
rates. A sinusoidal terahertz field applied to a semiconductor laser heats 
the electron-hole plasma and consequently modifies the optical susceptibil- 
ity. This allows an almost linear modulation of the output power of the 
semiconductor laser and leads to a faithful reproduction of the terahertz- 
field waveform in the emitted laser intensity. 
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Semiconductor laser modulation and switching via the pumping current 
at GHz rates is based on the so-called relaxation oscillationfl. 2]. Physi- 
cally. relaxation oscillation describes the periodic energy exchange between 
the laser field and semiconductor, which is an interband process dominated 
by stimulated electron-hole pair recombinations at high pumping level and 
occurs on a nanosecond time scale. For many applications, modulation and 
switching at higher rates is highly desirable. While improvements in modula- 
tion speed due to quantum well (QW) structure designs and other means have 
been quite impressive, further improvements beyond tens of GHz seem to be 
more difficult. The fundamental limit is placed bv the slow interband carrier 
dynamics. While current modulation has the unique property that allows 
direct and simple modulation, other means to modulate and switch a diode 
laser at higher rates are under intensive investigation. Notable examples are 
switching due to non-equilibrium carrier dynamics under femtosecond optical 
excitat,ion[3, 4. 5]. 

In this letter, we report, an alternative scheme for modulating and switch- 
ing semiconductor lasers by employing THz-frequency electric fields. Recent 
years have seen a surge of interest on the interaction of TITz waves with 
QW's. and TIIz field generation by semiconductor heterostructures. More- 
over, our research is motivated by an experimental demonstration of heating 
of an electron gas by applying a TIIz driving field (henceforth denoted by 
TDF) [6]. The experiment demonstrated that even under strong TIIz irra- 
diation. the electron gas maintains essentially a heated, thennalized form. 
Subsequent theoretical works [7. 8. 9] showed that the heated plasma has an 
almost sinusoidal temperature variation with time with a frequency given by 
twice that of the TDF. Additionally, the associated gain and refractive index 
change due to the plasma heating effect induced by a half-cycle TIIz field can 
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also lead to efficient transient optical modulation and may form the basis for 
an ultrahigh-speed all-optical switch [9]. 

The purpose of this work is to demonstrate that a continuous TDF can 
induce a modulated power output from a semiconductor laser with a temporal 
modulation that faithfully follows the driving field. We consider a standard 
edge-emitting semiconductor laser with a QW active medium. In addition to 
a standard QW edge-emitting laser geometry, a TDF is applied such that the 
polarization is in the plane of QW layer. Such a TDF can be applied from 
a free-spaee THz source or by integrating the laser with a photoconductor 
and THz transmission structure. Our theoretical approach starts with the 
following equation for the carrier distribution function 
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where a = e(+k),/?(-k) for electrons and holes, respectively, and F T h z is 
the TIIz field. The first two terms in the right hand side of Eq. (1) repre- 
sent carrier-carrier (c-c) and carrier-phonon (c-ph) scattering, while the third 


term accounts for stimulated recombination and the last term for pumping 
and other decay processes. We note that the TDF creates highly anisotropic 
carrier distributions (n* s 7^ r>k y ) which must be solved for on a two dimen- 
sional mesh of k points. Though a self-consistent many-body calculation of 
such anisotropic scattering terms is feasible, it involves a significant compu- 
tational effort, and thus we restrict ourselves to the standard relaxation-rate 


approximations for both of the c-c and c-ph scattering in this paper. The 
relaxation times (inverse of the relaxation rates) are 100 fs for c-c and 600 fs 
for c-ph scattering. Many-body interactions are not considered in this paper 
and the related effects will be discussed in a more detailed investigation. For 


3 



this paper, the THz waveform is assumed to be sinusoidal. The results, how- 
ever, carry over to more general temporal waveforms that may be of interest 
for switching and modulation with a signal field. 

For the description of the semiconductor laser, we employ a well-tested 
rate-equation model except we add additional terms for the TDF-induced 
plasma-heating effects. Under the rate-equation approximation, which es- 
sentially amounts to ignoring the polarization dynamics of the medium, the 
basic equations take the form: 
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where the symbols have their usual meanings: J, w, and P are the pumping 
current, active medium thickness, and the confinement factor; K, n g , and tf, 
are the optical wavevector, group index, and background dielectric constant; 
k is the cavity relaxation constant; E is the slowly varying amplitude of the 
laser field, and N = 2 Ylk n t/V * s ^ ie carrier density for electrons and holes 
which are assumed the same. The optical suscept'bilit.y, \ = \' + i\" , can 
be written as (linear approximation) [10]: 


vOV, T) = x(N (h To) + - ;V 0 ) + ^(T - T 0 ) . (4) 

We assume that the effects of the constant \'(/V 0 , T 0 ) are absorbed into the 
cavity decay constant k. The above expression can subsequently be cast into 
the following form by introducing the density and temperature a-factors [10]: 

\(iV. T) = -uiy(N - A7 () )(l - in N ) - in r (T - 7 0 )(1 - ia T ) . (5) 
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A few words about the model are in order. To start from a more fun- 
damental level, we need to solve the Boltzmann equations coupled to the 
optical field and TDF as well as to the optical polarization. This would cor- 
respond to the semiconductor Bloch equations (SBEs)[ll] with an additional 
THz field. The coupling of such SBEs with the wave equation for the optical 
field will offer a more complete framework for studying the interplay between 
the TDF, optical field, charge carriers, and phonons in a semiconductor. For 
lasers, a systematic approximation along the line of moment equations would 
be more feasible. Such an investigation is currently being pursued and the 
results will be reported elsewhere. 

Our present approach, represented mainly by Eqs. (2)-(3) can be justified 
in the following way. First, we ignore the heating effect due to stimulated 
emission of optical photons, as this heating is much less important (or lasers 
slightly above threshold than heating due to the TDF. This allows a sim- 
plification of the coupling between the heating and lasing processes; heating 
due to the THz field will influence lasing, but not visa versa. Under this 
approximation, the set of moment equations will be reduced to Eqs. (2) and 
(3) if we further assume that the relaxation of the average momentum and 
energy is much faster than the relaxation of the total carrier density and 
laser field. Due to this simplification, we can solve Eq. (1) for a given carrier 
density separately. The resulting distribution will be fitted to a Fermi func- 
tion which thus determines the plasma temperature at a given density under 
the influence of the TIIz field. Such a calculation allows us to obtain the 
temporal variation of the plasma temperature on the time scale of the TDF. 
The underlying assumption is that carriers after accelerating will thermalize 
quickly through c-c and c-ph scattering and a well defined temperature will 
oscillate with the TDF. This is a valid assumption due to the subpicoseconcl 



scattering processes and is also established experimentally in Ref. [6]. 

Figure 1 shows typical examples of the temperature evolution by appli- 
cation of TDFs with an amplitude of 2.5 kV/cm (circles) and 5.0 kV/cm 
(dashed lines), respectively. We see that after a rapid initial temperature 
increase from 30 I< (our chosen lattice temperature), the plasma tempera- 
ture is stabilized within ~1 ps to a sinusoidal form. The time dependence 
of the temperature as a solution of the Boltzmann equation after the initial 
transients can be fitted to the form: 


T = T 0 (F r// .) + T am (F THz ) cos(4ins Q t + <p 0 ) (6) 

where u Q is the frequency of the TDF. As shown by the solid line in Fig. 1. 
the fitting is almost indistinguishable from the actual data. Similar behavior 
for a stronger field (F 0 = 5.0 kV/cm) is depicted by the dashed curve. This 
time dependence of the plasma temperature is consistent with those seen in 
another theoretical work [8]. The T 0 and T am values are similar to those 
measured experimentally and calculated theoretically. 

In the following, we assume that the time dependence of the plasma 
temperature is given by Eq. (G). To study the effects of the time-dependent 
plasma temperature, we solve Eqs. (2)-(3) numerically for a typical edge 
emitting diode laser. As a first example, we consider a laser of 200 gm in 
length with cleaved facets at both ends. The power reflectivity is assumed 
to be 35 % for both facets. The active region consists of 3 QW's of 10 nm in 
width. The values of the a-parameters are taken from Ref. [10]. In all our 
examples, we use the values of T 0 =130 K and T am — 40 K. Figure 2 shows 
the laser intensity (solid lines) and carrier density (dashed lines) for two 
cases. In Fig. 2(a), F THi! is modulated at 50 GII/. while the corresponding 
frequency of the laser intensity is 100 GHz. We recognize from the figure 
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that there is a regular sinusoidal oscillation with a deep modulation, while 
the carrier density shows only a slight wiggle. In Fig. 2(b). the frequency is 
increased to 0.5 THz. The intensity yields the same regular oscillation, but 
the modulation depth is significantly smaller than in Fig. 2(a); additionally, 
the carrier density is now a constant over many oscillation cycles, indicating 
that the slow interband recombination is not involved in this fast modulation. 

With a further increase in the frequency of the TDF, the same regular in- 
tensity oscillations are observed. The modulation depth, however, decreases. 
For a frequency of a few THz, our model needs to be modified and will be 
discussed in a future investigation. In Fig. 3 two examples of low frequency 
TDFs are shown. As the temperature oscillation frequency decreases, the 
coupling with the carrier dynamics becomes important since there is a near 
resonant interaction between the plasma temperature oscillation and the re- 
laxation oscillation which is of the order of a few GHz. Figure 3(a) clearly 
shows signatures of such an effect. Consequently, the intensity-oscillation 
pattern is modified and the carrier density changes appreciably during the in- 
tensity oscillation. For the slightly higher temperature-oscillation frequency 
in Fig. 3(b) (25 GHz), we see a deeper modulation. Now the laser is almost 
switched off. The intensity pattern of Fig. 3(b) also shows a period-doubling 
phenomenon. This is again caused by an interplay between the carrier and 
temperature dynamics. 

To conclude, we have illustrated an ultrafast scheme to modulate and 
switch semiconductor lasers at THz and sub-TIIz frequency rates. This is 
possible by circumventing the usual slow interband carrier dynamics. We 
show that under a sinusoidal TDF. the temporal pattern of the electrical field 
can be converted to a similar, quite regular pattern of the laser intensity. This 
faithful representation of the sinusoidal THz field by a laser intensity pattern 



in the visible or infrared domain may offer an alternative way of measuring 
high-frequency alternating electrical fields. For lower-frequency oscillation, 
the same arrangement can lead to an almost 100-% modulation and even 
on-off switching behavior. This could be useful for switching semiconductor 
lasers. Currently we are conducting a more detailed study using the moment 
equations mentioned above. 
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Figure captions 

Fig. 1 Plasma temperature oscillation hv means of a terahertz driving field at 
two field amplitude values: 2.5 kV/cm (circles) and 5.0 kV/cm (dashed lines). 

Fig. 2 Laser intensity and carrier density versus time in the presence of a 
terahertz driving field of frequency 50 GHz (a) and 500 GHz (b). 


Fig. 3 As in Fig. 2 but with a field frequency of 5 GHz (a) and 25 GHz 


(b). 
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